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A SPECIAL METHOD FOR FINDING BODY DISTORTIONS THAT REDUCE THE WAVE DRAG 
OF WING AND BODY COMBINATIONS AT SUPERSONIC SPEEDS 1 


By Harvard Lomax and Max. A. Heaslet 


SUMMARY 

For a given wing and supersonic Mach number, the problem 
of shaping an adjoining fuselage so that the combination will 
have a low wave drag is considered. Only fuselages that can 
be simulated by singularities ( multipoles ) distributed along the 
body axis are studied. However, the optimum variations of 
such singularities are completely specified in terms of the given 
wing geometry. An application is made to an elliptic wing 
having a biconvex section, a thickness-chord ratio egual to 0.05 
at the root, and an aspect ratio equal to S. A comparison of 
the theoretical results with a wind-tunnel experiment is also 
presented. 

INTRODUCTION 

Tlie most simplifying assumptions that still permit the 
construction of a mathematical model general enough to 
contain quantitative information about steady three- 
dimensional supersonic flow are those basic to the develop- 
ment of linearized theory. Of these, the two principal 
assumptions are that the viscosity effects are negligible and 
the perturbation velocities are almost everywhere small 
enough to be neglected relative to the flight or free-stream 
velocity. Under such restrictions the flow field can be 
described in terms of a perturbation velocity potential <p 
obeying the equation 

P 2 <Pxx—<Pn—<Pu=0 (1) 

where 1 and the reference coordinate system 2 is 

shown in figure 1. Further, the wave drag of any object in 
a flow field governed by equation (1) can be evaluated (see, 
e. g., ref. 1) by means of the equation 

D= — po f dx f dO\lim (2) 

J ~eo J 0 L r— _1 

where x, r, and 6 are cylindrical coordinates also defined in 
figure 1. 

General solutions to equation (1) are numerous and clas- 
sical. In applying these solutions to the interpretation of 
physical phenomena the usual approach is to fit them to the 
given boundary conditions, that is, to make the flow field 
simulated by them conform to the shape of the disturbing 
object as well as to a uniform free stream at infinity. Hence, 

i Supersedes NAOA RM A66B16 by Harvard Lomax and Max. A. Heaslet, 1065. 

* It should be stressed that the x axis is parallel to the free-stream direction (wind axes) 
ao a body of revolution can be symmetrical about this axis only at tero angle of attack. 


from this point of view, the choice of a type of general solu- 
tion to be used in analyzing a particular problem with the 
least mathematical effort depends on the geometric form of 
the object under consideration. For example, general solu- 
tions based on Green’s theorem are well adapted to the study 
of forces on single planar wings in a steady supersonic flow. 
On the other hand, the general solution given by Lamb (ref. 
2) — which is composed of an infini te set of multipole distri- 
butions disposed along a line — is well adapted to the study 
of the flow around fuselage-like objects. 

In this report use is made of certain general solutions to 
equation (1) but with a deviation from the usual approach 
mentioned above. One considers, in fact, two different 
kinds of solutions which represent separately, in a given 



Figure 1. — Reference coordinate systems. 
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vicinity, different classes of real objects and, by means of 
equation (2), finds optimum combinations of these solutions 
from the viewpoint of low wave drag. The analysis involved 
in solving this problem has, in general, a distinct mathemati- 
cal advantage over the problem of calculating the drag of a 
given object; namely, that the immediate problem of finding 
a shape with a relatively low wave drag is divorced from 
any detailed reference to the shape itself. It is true, of 
course, that the stream surface representing this shape must 
eventually be found and, in fact, a limi tation on the appli- 
cability of the method is given by the requirement that this 
shape be real. However, the problem of finding the shape 
of the object when <p is known is a matter of direct calcula- 
tion. 3 One should also be careful to notice that the op timum 
solutions obtained by this procedure are not necessarily true 
optimums but purely relative to the choice of solutions used 
in the analysis. 

LIST OF IMPORTANT SYMBOLS 


A 

A n (x) 


a„(x) 

B n (x) 

b 

».(*) 

C D 

c, 

D 

D n 


M 

<J 


r t 

R 

S 

S o (x,0) 

t 

TJo 

V 

x,y,z 

x,r,e 


wing aspect ratio 

( — 0) n times the nth derivative of the nth multi- 
pole distribution a n (x) (See eq. 16.) 
semi-root-chord of elliptic wing 
strength of nth-order multipole distribution 
multiplying cos nO 

(-— /3) n times the nth derivative of the nth multi- 
pole distribution b K (x) 
semispan of elliptic wing 

strength of nth-order multipole distribution 
multiplying sin nO 

drag coefficient, 

pressure coefficient, local pressure minus static 
pressure divided by g_ 
wave drag 

wave drag associated with nth-order cancella- 
tion multipole distribution (See eq. (59).) 
maximum fore-and-aft extent of wing equiva- 
lent multipole distribution 
maximum fore-and-aft extent of wing equiva- 
lent multipole distribution for angle 6 
free-stream Mach number 

free-stream dynamic pressure, p ° - 

m 

L„ 

see equation (46) 
radius of body 
area of wing plan form 

normal projection of wing cross-sectional area 
measured in oblique planes 
maximum thickness of wing root chord 
speed of free stream 
volume 

Cartesian coordinate system, x parallel to free- 
stream direction 

cylindrical coordinate system, x parallel to free- 
stream direction 


* From a mathematical point of view the essence of the method outlined above Is that the 
analysis Involves the solution to direct problems, that is, problems of Integration. Cal- 
culating the drag of a given body, on the other hand, involves the solution to inverse prob- 
lems, that is, problems involving the inversion of integral equations. 


a n (x ) strength of nth-order optimum cancellation 

multipoles (See eq. (36).) 

. /3 

X„ slope of wing upper surface .measured relative 

to free-stream direction 
n tan -1 (/8 cos 6) 

p 0 free-stream density 

ip perturbation velocity potential 

DEFINITION OF THE PROBLEM 

The problem of designing an airplane to have a minimum 
wave drag must be stated quite precisely. If the aerody- 
namicist is approached with the question, “Given an aero- 
dynamic shape, can its wave drag be lowered?” he can always 
reply that any volume of material having a wave drag can 
always be reshaped within a space of finite dimensions so 
that it will have less wave drag at a given Mach number. 
Such an answer is interesting but, at present, not very useful 
to the airplane designer. There is first, of course, the basic 
criterion that the total drag should be minimiz ed at a given 
. lift and minimizing a component part of this total without 
holding the other parts fixed does not necessarily yield the 
lowest possible drag for a given set of restraints. For ex- 
ample, the configuration illustrated in figure 2 has no wave 


Cylindrical shroud 



Mach waves 

Figure 2. — Body and shroud with zero wave drag. 


drag when traveling at zero angle of attack; but it has a 
relatively high friction drag, because of the largo amount of 
wetted area, and its drag due to lift could also bo relatively 
high. Completely aside from all such performance consid- 
erations, however, are many other important considerations 
that are unfortunately more or less vaguely defined from an 
aerodynamic point of view. For example, an airplane must 
contain a certain amount of usable volume, the shaping of 
individual parts is limited by structural requirements, and 
the arrangement of these parts must not seriously harm the 
airplane stability and control. The interrelation of all such 
separate demands presents an extremely complex design 
problem making it difficult to deviate too far from the reliable 
shapes set by experience. 

As a result of the above-mentioned difficulties, the aero- 
dynamicist who is concerned with discovering a practical 
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airplane shape having low wave drag finds the real definition 
of his problem somewhat obscure. In a sense his first prob- 
lem is, literally, to pose a problem; that is, to impose a 
minimum number of arbitrary but pertinent restraints within 
the framework of which the wave drag is to be minimiz ed. 
Even when this has been done, he still is concerned with the 
question of uniqueness, since optimum shapes are not neces- 
sarily unique even when several restraints are imposed. 
Consider, for example, the problem of finding the Busemann 
biplane which will have minimum wave drag at a given 
Mach number for a fixed section strength, volume, and 
wetted area. If the design Mach number is 1.41, one such 
design (on the basis of linearized theory) is shown in figure 
3 where the chord-gap ratio, h/c, is equal to %. The resulting 
variation of the wave drag is shown in the upper part of 
figure 4. However, when the gap is closed to the point 
where h/c equals %, the variation of wave drag, shown in 
the lower part of figure 4, is the same within the interval 
1.28 <M< 1.66 and everywhere else is lower. It is likely 
that one would have first discovered the former solution, yet 
to the accuracy of the theory used, the latter is obviously 
preferable. 

With the above observations always in mind, attention 
will be directed in this report to the analysis of simplified 
configurations composed of two distinct types of volume: 
planar types, that is, wing-like volumes, thin in one dimen- 
sion and bounded by surfaces that never deviate far from 
a reference plane; and rectilinear types, that is, fuselage-like 
volumes longer in one dimension than in the other two and 
disposed more or less symmetrically about a straight line. 

In particular the following problem is posed: 

Given a thin nonlifting wing, what is the shape of an 
adjoining fuselage, the stream surface of which is 
simulated by a line of multipoles in the same plane as 
the wing, that will minimize the wave drag of the 
combination at a given Mach number? 

BASIC CONCEPTS 

A LINE OF SOURCES 

The velocity potential induced at the point x,r,d by a group 
of sources distributed along the x axis, starting at —L 0 , is 
well known to be given by 



Figure 4. — Drag variation for double-wedge Busemann biplanes. 

where a„(£) is the source strength per unit of length. In 
order to calculate wave drag one needs only the value of v as r 
approaches infinity. This asymptotic value is simple enough 
to find provided it is observed that, as r is increased, x should 
also be increased so the potential can be studied in the 
vicinity of the Mach waves radiating from the disturbing 
object. Hence, set 

x=x 0 +p r (4) 

so for a given r, x 0 measures the streamwise distance of the 
point x,r,6 from the Mach wave emanating from the origin 
and, in particular, the foremost wave is located at x 0 = —L 0 . 
(See figure 5.) 


<p(x,r,0)= 


2 TV 


"t-er a »(&% 


(3) 


Mach waves at M = 1.41 



Figure 3. — Busemann biplane. 
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Figure 5. — Coordinates introduced in equations (4) and (6). 
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If equation (4) is placed into equation (3) and r is assumed 
to be large, the potential induced by a source is 


v)r 


_1 C x ° Qp (£)d£ 

-J28r J 4x6— 


2 ir-yj2pr J-io4x 0 — £ 
and the induced velocities are 


<f>i\ 


_i 

■\/28r J -A> Jx„— 


2 TrV^jSr J~ L o 4^o— k 

MULTIPOLES 


(5) 

(6a) 

(6b) 


Lamb, in reference 2, page 527, has presented a general 
solution to equation (1) consisting of an infinite set of basic 
singular solutions. These basic singularities, referred to as 
multipoles, can be distributed along a line and weighted so 
as to reproduce certain body shapes enclosing the line. The 
expression for the perturbation velocity potential for a 
distribution of 7ith-order (n= 0,1, 2,. . . ) multipoles starting 
at — L 0 and continuing along the x axis can be written in 
terms of a cylindrical coordinate system (fig. 1) as 


r'/l by [«*(£) cos 7i0+ &,(£) sin 7i0] 

^’ r ,0 )= ~2^Vr W J-l. “ 




(7) 


as 


The operator ^ is defined < 

(1 d_V = l il/l+ 1 d 1 d* 

\r brj r dr \r br) r 3 br r 3 dr 2 

/I d\* 

and the definition of ( - follows by induction. If the 
notation 4 

{AX f z A A) d v ( ^„ (i)(3) ... (2 ti-i) C x A(v)dy 
W J a 4 ^ K l) 2“ J a 

is introduced, where the symbol J- is read “finite part ot the 

integral,” equation (7) becomes 

/ rr M r"P’ Jn ( 2 n)\ jy-P' Kfe) cos 7i0+ 6,(£) sin 7i0]d£ . . 

<PnV,r,u) 2' l + 1 +»! f- L „ [(x-tf- 1 SV] (S ' +1)/J W 

and the general expressions for the induced velocities be- 
come — writing only the term involving the cosine, since the 
result for the sine is directly analogous 


JrrflX r"fl*»(2ft+2)! f*+e r (s-+)q,(£)cos7 iBdB 
^ 2 2 "+ Vii + l )! T~^ [( r -^- iSV 2 ]' 2 "- 3 ^ 2 


<Pr(X,r,d)=—J2 


P2n r n-l(2n)\ fx-flr 


^ 2 n+1 n(ri)l T-to 

[7t(x— SP+^+l^V]®^) cos7i0d£ 

[(z-S) 2 -/3 *r*] 

1 (rr rtr n ~ l ^ in (2n)l Jj-sr a,(£) sin7 i9d£ 

r m¥, J Ai 2 a+l Tr(n)\ J- fo [(x-$) 2 -iS+] (J «+ 1 >/ 2 


(9a) 

(9b) 

(9c) 


Another very useful way of developing these multipole 
solutions evolves from an application of operational tech- 

• For a detailed discussion of the finite-part concept as used in this report see reference 3. 


niques. To begin with, rewrite equation (1) in terms of a 
polar coordinate system, thus 


dy dy i ay 

p bx* br 3 rbr r 3 be 3 
Next, define the Laplace transform of <p(x,r,d) by 

r* oo 

£(s,t-,0)= I <p(x,r,e)e-"dx 


(10) 


(ID 


and apply this transform to equation (10). There results 
(for a proof see Appendix A). 


„22-_ay_i dy 

P v br 3 r br r 3 be 3 


( 12 ) 


Now, if a general solution to equation (12) is expressed in 
the form 

v(8,r,d)=f(r) cos nd 
then /(r) must satisfy the equation 

Solutions to this are given by 

f(r)= 1 7 n {8)Wra)-\-A n {8)KJfir8 ) 

where J, and K n are modified Bessel functions as defined in 
reference 4, page 77. Hence, if <5 is to vanish as r goes to 
infinity, a general solution to equation (10) can be written 
in the form 

7p{s,r,e)=~ {A, [X,(s) cos 720+I?„(s) sin n9]K n {fir8) (13) 

67T 0 

The above result will be transformed back to the physical 
plane in two ways. First, apply the identity (ref. 4, p. 79) 




and re-express equation (13) as (only the coefficient of the 
cos nd term is written since the treatment of the sine term is 
identical) 

^ r , 0 ) = __L s (- 1 ). cos ne (} £)* Kotfrs) 

The inverse Laplace transform of K 0 (firs ) — see reference 
5 — is 

0 


L~'[K 0 (firs)\= 




x</3r 

x>/3r 


So, since 

j_ L dx? . . -Jit' ^«A n (x„)= 

( 5 = 1)1 A. , ( 14 ) 
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an application of the convolution integral and other standard 
operational techniques yields 


V( x-tf-pv 


i \rdrJJ-i V(*-f)*-|S*r* 


: 


(15) 


From comparison of equations (8) and (15), the relation 
between the strengths a n (x) and A n (x) for the two different 
forms of the solution is found, to be 


n=0 


a 0 (x)=A 0 (x) 

(— /3) s a,,(z )=— L (*— xJ'-'AJxddxi, n>0 


or 


(-pya?(x)=A n (x) 


(16a) 


(16b) 


where (x) symbohzes the operation a H (x) and where 
use is made of the conditions 

(-£.) =tf- u (-Z 0 )= . . . =a n (-L o )=0 (17) 

Another way to transform equation (13) back to the 
physical plane is to do so directly. In this way one finds 
(from ref. 5) 


L 1 K n (firs)= 


L^ttl 


cosh | ft cosh -1 ( — 
ir 


, x<P r 


x>pr 


from which equation (13) reduces immediately to 


<p(x,r,6)= — jj- X cos nd 


X 


t-prA^) cosh Ticosh - ^^^) 


£o 


V(i-f) 2 -/9V 


(18) 


The perturbation velocities in the field represented by 
this potential are readily calculated. Thus 


<p x (x,r,d)= — 7 — X) cos n0 


1 1 " 

- ^(r,r,6>)=— X) » sin nd 

r Zirr 0 


J r-0r 
-£0 


A,'(Q cosh n cosh -1 (^r)J 


V(x-?)*-i 9 *r* 


x; 


(19a) 

■*-0, A ({) COSh COSh -1 (^r)] 


V(r-Q 2 -/3V 


(19b) 


and by taking the derivative of equation (13) with respect 
to r, one finds 

psA o {s)K0rs) +| X fisA n (s)[A:,_ 1 Ors)+ K n+1 (prs)] 
which transforms to 


r” ■- r- ^- ,(a {^ V cosh[»coA-’(y)]}« \ 

ir\J-Lo r j(x-Zy-P 2 A 2 i J-Lo y(*-0*-0V / 


(19c) 


If the relation between the functions a n (x ) and A n { x) is 
given by equation (16), the velocities represented by equa- 
tions (9) and (19) are, of course, identical. 

In order to obtain limiting values induced by multipoles 
distributed along the x axis starting at —L 0 , one returns to 
either equation (8) or (18) and calculates the leading term in 
a 1/r expansion. As in the derivation of equation (5), it is 
necessary to observe that as r is increased, x should also be 
increased so ^ is given in the vicinity of the foremost Mach 
cone created by the multipole distributions. Hence, using 
equation (4), one finds for equation (18) 


<p=~— X) cos 710 

4TV 0 


X 


A,($) cosh 

T 1 

0 

0 

g- 

1 

iff \ x ° 


V 11 Pr ) 

' V 

2/3r(x 0 -Z)t 

l 1+ *) 



which has the leading term as r goes to infinity 

\ 1 ^ „ p° A n {m 

<p) r-.m = 7=- > , COS nd 1 

2^^^f2f}ro J-x 0 jx 0 — £ 


( 20 ) 


Similarly, the perturbation velocities reduce to 


#r)r— ►<» == 7= 2D 008 716 

2 Tnj2fir 0 

^r)r— ►<»— fi<p x 


X 


*• A n \m 

a-? 


r- Ajm 

J~ L o A>— £ 


(21a) 

(21b) 

(21c) 


In calculating the wave dr-ag using equation (2) only the 
velocity components <p x ) r - » „ and ^ r ) r _ m are necessary. 
Hence, from comparison of equations (21a) and (21b) with 
(6a) and (6b), it follows that at a given 6 a series of multipoles 
induce the same momentum flux on an infinite cylindrical 
control surface as a line of sources having a strength varia- 

CD 

tion a 0 {x) equal to X cos nOA H '(Q. If one identifies a line 
0 

of sources with a body of revolution, then it is apparent 
that, at a given 6, a dragwise equivalence has been estab- 
lished between a line of multipoles and a body of revolution. 

HAYES’ THEOBEM AND ITS APPLICATION 

In the previous section a relation was found between 
multipole and source strengths which produce, at a fixed d, 
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z, 



Figure 6. — Mach forecone from x,r,B in xi,yi,zi space. 


equivalent momentum transport across a cylinder of infinite 
radius. By using a theorem attributable to Hayes (ref. 6 ) 
one can derive the strength relationship between any distri- 
bution of singularities throughout space and a line of sources 
which gives the same equivalence. 

The essence of Hayes’ theorem is that, for a fixed B, the 
velocities induced on a cylinder of infinite radius by singular 
solutions to equation ( 1 ) (e. g., sources and doublets) are 
invariant to displacements of the singularities along certain 
oblique planes. In order to be specific, the equation of 
these oblique planes is next derived. 

Consider the point x,r,B in a flow field having a supersonic 
free stream moving parallel to the x axis. Figure 6 shows 
the Mach forecone (by definition the Mach forecone is the 
boundary of the region within which a disturbance in a 
supersonic stream can affect the flow at the cone apex) 
from x,r,B in x^y^Zi space. The equation of the forecone is 

Xy=x~p V(r cos 0 — 2 /i) 2 +(r sin B—ztf (22) 

One 'wishes to let r become very large and find the shape of 
the forecone as it passes through regions close to the origin 
of the Xi,yi,Zi coordinate system, regions in which the objects 
creating the wave drag are located. From equation (4) 
and the expansion of equation ( 22 ) for large r, it follows that 

Zi=£o+/3r — pr |\-~ (j/i cos d+z t sin 9) + Vl 

£} 

=Xo+/3(yi cos B+Zi sin 6 )— ^ ( y x sin 6— z l cos fl) 4 -!-. . . 


and when r goes to infinity 

£i=£o+/3yi cos 6+fizi sin 6 (23) 

which is the equation of the oblique plane mentioned above. 
It should be noted that the envelope formed from these 
planes by fixing x 0 and varying 6 between 0 and 2 -jt coincides 


with the Mach forecone and aftercone from the point x 0 , 0,0. 
Hayes’ result can now be stated 6 

To the lowest order in 1/r, as r tends to infinity, the 
magnitude of the perturbation velocity potential and 
its gradients at a fixed angle 6 is invariant to a finite 
translation of sources (or any other singular solution 
to the wave equation) on planes parallel to that given 
by equation (23). 

Consider the application of Hayes’ theorem to planar 
distributions of sources lying in the Zi=0 plane. As is well 
known, such a distribution simulates a wing symmetrically 
disposed about the horizontal (zi=0) plane. In fact, if 
X«( 2 i,yi) is the local slope of the wing upper surface, the 
local source strength per unit area (according to thin airfoil 
theory) required to simulate the wing is — C/ 0 X m / 7 t and the 
velocity potential of the disturbed flow field is given by 


U» P C K(^i,yi)dx l dy 
* V(x— Zi) s — /3 4 (y— yi)*— 0 a z 2 


(24) 


where r is the area of integration bounded by the wing edge 
and the trace in the zj=0 plane of the Mach forecone from 
the point x, y, z. Next introduce the new coordinates and 
yi such that £i lies along the Xi axis and ??i lies along the 
intersection of the Zi =0 plane and the plane given by 
equation (23) (see fig. 7). Set 


and 


tan a=/3 cos 0 
£i=*r- 2 /i tan /tl 


(26) 


sec n 

Zi=£i+ 77 i sin n 


(26) 


2/1=771 COS n 


Then, in terms of the $ 1,771 system, equation (24) becomes 


Uo f f ^n(£i+??i sin y ,, 771 cos n) cos n dy 1 

* J'J V(^ — — ^71 sin cos nY-Fz 2 



« For proofs, see Hayes’ original derivation (ref. 0) or, if more convenient, reference 1. 
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Wing plan form-. 



Figure 8. — Position of wing sources and equivalent single source. 

As before, the asymptotic value of p as r=Vj/ 1 +2 2_> “ is to 
be calculated. Accordingly, one can apply Hayes’ theorem 
and sum up all the sources along a line £ 1 = constant (e. g., 
between a and b in fig. 8) and place them as a single source 
on the axis. The strength of this equivalent single source 
is'given by 

Z7 o <S' B ,'(£i,0)=2t7 o cos /i X„(^+ ij t sin /*,ui cos n) dm (28) 
J tcing 

where the integration is taken across the complete wing along 
the’line {^constant and S w '(^,0)=^ Sd£i,9). 

The term S u (x,9) has a clear geometrical interpretation 
(see fig. 9), being simply the normal projection of the wing 
area intercepted by the oblique plane 6 Xi=x-\-f3yi cos 9. 

The above defines the strength variation of a line of 
sources (and, therefore, a body of revolution) which induces, 
for large r and a fixed 6, a potential field identical to that 



Figure 9. — Wing area intercepted by oblique plane. 


• The true oblique plane la given by equation (23) but the wing is so close to the n«0 plane 
that the variation with zi can be neglected. 


induced by a given wing. Hence, the results given in 
equation (6) yield 


Up r*o i S' j'ikifi) 

2tt V2/3r J — 1 ,( 8 ) V x o ?i 


(29a) 


<p r )r->a> = fitpx (29b) 

A similar result exists for a planar doublet distribution 
(see ref. 1 or 6) but, in this report, only problems in which 
the wings have no loading (local lift) will be considered. 
Lifting effects have been treated in a similar fashion in 
reference 7. 


CANCELLATION MULTIPOLES AND DRAG MINIMIZATION 

Since the flow field is governed by a linear partial differ- 
ential equation the velocities induced by different solutions 
to it are additive. Therefore, the drag of an object simu- 
lated by various multipoles distributed along the sq axis 
and a sheet of sources in the Zi=0 plane is given by 

D = ~PoJ o J" dx^ lim r t((p r )„+(95 r ) J [(?>*)«,+ (p*)»] 

where the subscripts m and w refer to the multipoles and 
wing sources, respectively. But equations (21) and (29) 
identify, for a fixed 6, these' velocities with those induced by 
equivalent line sources. Hence, for any given 9; one can 
immediately apply K&rm&n’s drag formula (ref. 8) and 
then for the total drag, integrate 9 from 0 to 2v. This leads 
to 


„ f (‘L (8) PL (8) 

D=-£k \ dB \ dx ' 

OTU Jo J_£'(0) J-L'(9) 

^5' tt "(x I ,0)+S(-|3)»[a?+ 1 > Oncost W+b« + » (*,) sin nflj- 
{tWW)+^ (-mat +v (a*) COS n9+ 


b J +U (x 2 )sinn^] 


J-lnlaq— aql 


However, since both the wing and multipoles are in the same 
plane (interpreted physically, the wing is centrally mounted 
on the fuselage) and the wing is simulated by sources only 
(has no twist or camber), one can show the optimum value 
of each b ? +1) (aq) is identically zero. 7 Hence, one can write 


„ r 2r rL (8) fL (8) 

D=~B\ dB *q dx i \UoS 0 "(x u 9)+ 
OK* Jo J-Z'(6) J-Z'(0> L 

f) (— (3)*a. ( * +1) (aq) cos nil j 


U 0 S a "(Xi,9)+iti (— /3)*®? +1) (zs)cos ni\ ln|aq— aq| 

o J 


( 30 ) 


7 By symmetry 
Hence 


J. 


S w "(r,tf)-S."C r,2*-0) 

lr 


sin n3Sm"CrJf)dQ^0 


and any positive or negative variation of b n C«H) (n) can only increase the drag. 



716 


REPORT 1282 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Next expand the term S w "(x,6 ) in a Fourier series, 
finds 

One 


U 0 S„"(xJ)='y' l (—B)*a < ? +1) (x) cos nd 

Q 

(31) 

where 

£ T SJ(x,e)de 

(32) 


(*)=^z^ j o S w '(x,6) cos nBde 

(33) 


Place these expressions in equation (30), integrate with 
respect to & — using the orthogonal property of the trigo- 
nometric series — and one finds 


D=2D 0 J r'^2iD n (34) 

i 

where 


D„ 


/5 2 


P dx i p <M«/ n+1 H*i)+a. ( " +1) (zi)] 

o J-L„ J -L„ 


g iirU^ 

[«. < " +u (^)+a, ( » +u (* a )]l7ijx I -%| 

Since one can show 8 

/=— f dx if dxj'{xi)f' {x^ln^—x^ >0 
J -L„ J 


(35). 


the minimum value of D as expressed by equation (34) is 
given when each D n is itself a minimum. In other words, 
each D n can be minimized separately. Further, it follows 
that the value of the minimum itself is zero and occurs when 


a ^+»(x)=-a n ^+»(x) 

f TT rtr 

-~Jo s u "(x,6)de, , n =0 

H rj rir (36) 

, - -s, S (x,d) cos n0dB, n > 0 

s. p) J 0 

Equation (36) is the mathematical definition of the 
optimum cancellation multipoles; namely, those multipoles 
which are just equal in magnitude and opposite in sign to 
the wing equivalent multipoles — equivalent in the sense 
that they induce an identical momentum flux across a 
cylinder of infinite radius. 

Obviously, if all the optimum cancellation multipoles 
were used, the wave drag of the combination would be zero. 
This result must, however, be properly interpreted with 
regard to the simulated shape. In order that the multipole 
lines can represent the distortion of a real fuselage, one 
must assume a cylindrical body exists upstream from the 
Mach cone x-\-L 0 —fir (the effects of the nose are being 
neglected). This body forms the initial boundary of the 
stream tube which represents the physical fuselage in the 
vicinity of the wing and multipole lines. Clearly, the area 
enclosed by this initial boundary can be small enough for 
the subsequent stream surface to cross itself and represent, 
therefore, a physically unreal body. Hence, the fact that 
the wave drag of the wing and multipole combination can 

00 j CD 

• Sct/Ci) = 2 A. sin rut, x= —L. cos Integration gives I=^r 2 "AV which can never 

1 1 1 

be negative. 


be reduced to zero is quite valid, but in the over-all picture 
not only have the inevitable nose and tail drags been neg- 
lected but also the shape simulated by the combination can 
be unrealistic. 

SOME PROPERTIES OF THE CANCELLATION MULTIPOLES 

Let us consider next some of the restrictions necessarily 
imposed on cancellation-multipole distributions and some of 
the particular properties of those given by equation (36). In 
the first place, if a M (x) is any multipole distribution that 
generates a potential field given by equation (8) or (15), it 
follows from equation (17) that the value of a,(x) and its 
first n derivatives should be everywhere continuous. Further, 
if a K (x ) is a constant behind some point, say L 0 (i. e., for 
°=)>x)>L u ), the induced flow field would simulate expanding 
streamlines in the case 7i=0 or some form of vorticity in the 
case 7i>0; the former case is to be avoided since any simu- 
lated body is assumed to have a finite area at x— <», and the 
latter case is to be avoided if there are no resultant forces 
normal to the free stream. 

One can show that all the above properties are satisfied 
by «,(£), the optimum cancellation-multipole distribution 
as defined by equation (36). First, notice that aj" +1) ({) 
must be zero everywhere outside the wing-enclosing Mach 
'forecone and aftercone, that is, for — co <£<— i/ and 
Z 0 <£< co (see fig. 10). (Any multipoles in these regions can- 
not combine with the wing equivalent multipoles and must, 
therefore, increase the drag.) Hence, one can set 

(— X 0 ') =«<— » (— i 0 ') 


= . . . «.(-Z»0 


=0 



Figure 10. — Symbols used in study of multipole properties. 
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Then tlio condition of continuity is automatically satisfied for 
ai B) (£) in the entire interval — « <^x<C . m if SJ(£,d) is 
derived from, a wing having finite wave drag (in particular, 
from a wing having no blunt edges along which the normal 
component of the free-stream Mach number is unity or 
greater). It follows immediately that a™ (£) , m<Cn, is 
continuous since the latter is found by integrating (further 
smoothing) (£) . 

The proof that a m (L 0 ) as given by equation (36) is zero 
requires more consideration. One can show, however, that 
“n m) (7/fl)=0 where 0 <m<n. First, the equality (£«,) =0 
follows from the fact that the wing closes and S a ' (L 0 ,9) itself 
is zero. Next consider the definition of Thus 

c03 7Wde 

— 2,u r r 

== >- ■ cos fi cos nd dd I X«(£+»? sin y,r) cos y)di] ■ 

( — P) 0 J A 0 (9,f) 


Further, this connection was direct and relatively simple if 
the strengths and positions of the distributions were given. 
Unfortunately the connection between the multipoles and 
the shape of the simulated surface is generally not so simple. 
Such a relation does certainly exist, however, and if the 
strengths of the multipoles are known, the relationship is 
again direct. That is, a given distribution of multipoles 
yields directly, by the formulas given in the previous section, 
the induced velocities everywhere in the flow field, these, in 
turn, fix the stream surfaces along any one of which (since, of 
course, the theory neglects viscosity) a physical surface can 
be imagined. 

In general, if 

F c (x,y,z )= 0 

F p (x,r,e)=0 

are the equations of a stream surface in cartesian and polar 
coordinates, respectively, then the equations 



where use has been made of equation (28) for the definition 
of S a '(£,0) and h 0 and hi are defined in figure 10 . Since 

, (?-? 1 ) n - m - , «n 00 fe) (37) 

one has 

(£)=; — o\n W ° i f 008 & 

(— /3)“iiT(m— n)J 0 


(tj I ZF C bF c 

( U o+<M d^ =0 

(tj , ^F ZF 1 ^ 

( Z7 ‘>™-5T+^"5T+^ v, V =0 


(40) 


must hold. 

For example, in studies of thin wings lying in a plane, the 
particular form of equation (39) 


J t r* i(«.ti) 

d,T}(t—%iY~ m ‘X^+ 7 / sin y,v cos y) cos y 

-V */*>(«, t,) 

Change the £ 1,17 coordinates back to the x,y system by means 
of equation (26) and this becomes 


z~h(x,y )= 0 

is assumed and equation (40) becomes 

— (Uo+Vx) a5~^^+^ = 0 


“ 3 da !j {i ~ x+fy oos e} — 1 

\ u (x,y)dxdy 


or, neglecting second-order effects, 


£>A_ 1 
bx~U 0 v ‘ 


The area S\, shown in figure 10 , becomes independent of 9 
when £=£„ (being then just the area of the wing itself), 
therefore 

a " W ( L °) = (-0)v r (m-n) jj^( x ,y)dxdy 
2 s 

J (L 0 — x-YPy coBd) n ~ M ~ l cosnB dB=0 


since, for m<yi 


which is the familiar boundary condition used in thin-airfoil 
theory. On the other hand, if the equation of the body 
shape is written in the form 

r—R(x,8 )= 0 

then equation (40) becomes, for linearized theory, 


. i)R <pe £)R 

0 ~bx~ <p ~W t 


(41) 



cos md coanB dd=0 


If the flow field is radially symmetrical or if the body surface 
is quasi-cylindrioal, equation (41) reduces to 


Hence, for the «*(£) defined by equation (36) 
a‘"> (L 0 ) = «*-” (L 0 ) = (Z„) = ... =a n (L 0 )= 0 (38) 

AIRPLANE SHAPE 

In the previous section a connection was established be- ; 
tween multipole distributions and their resulting wave drag. 


5 R 1 


(42) 


which is the familiar boundary condition used in the study 
of quasi-cylindrical bodies or bodies of revolution. 

In general, a nonlinear partial differential equation of the 
first order such as equation (41) can be reduced to two simul- 
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taneous, ordinary, nonlinear differential equations of the 
first degree (see, e. g., ref. 9). Thus equation (41) becomes 


dO 1 , r> 

dTuji^ x ' R ’ e) 

dx^Uo <pr( - x ’ R ’ e '> 



and if epe and <p, are known functions of x, R, and 8, these can 
be solved numerically. 

If the strengths of all rectilinear multipoles and source 
sheets are given, equation (19) or (24) can be used to find 
<p r and <pe at arbitrary field points. Hence, the first step in 
finding the body shapes reduces to that of integrating such 
equations. However, these integrations are difficult and 
tedious even when entirely numerical procedures are em- 
ployed and the results still have to be interpreted in terms of 
the body shape according to equations (43). Therefore, 
from a practical viewpoint, it is necessary to study certain 
approximate methods for obtaining the velocity field. 

Let attention be concentrated on the disturbances created 
by a line of multipoles. In particular, consider the fields 





_£• Exocl linearized theory 

t= ^ Slender-body theory 

00 Large 7 theory 



x o 


(b) Velocities, <p 0r . 
Figure 11. — Continued. 


induced by simple polynomial distributions satisfying, in 
each case, the end conditions given by equation (17). For 
particular variations set 

OoQc)=c 0 ( 1— a?)x ' 

a i (5)=c 2 (l— i*) 3 ^ - (44) 

a 4 (x)=c 4 (l— P)^, 

where c 0 , c?, and c 4 are constants determining the amplitudes 
and x=x/L 0 . Figure 11 (a) shows the variations of these 
coefficients with x, and figures 11(b) through 11(f) show 










2v<j> r • 2 tti j> r 

£ 3 e 2 cos 20 / . jS%- 2 cos20 
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(e) Velocities, <p tr . 
Figure 11. — Continued. 


how velocities induced by these distributions vary with 
%o (*«=* — /3 r) and r Cr=f3r/L 0 ). The results have been com- 
pared with those for large F given by equation (21) and with 
those for small F given by slender-body theory. Values for 
the latter theory are deter min ed from equations (9) or (19) 
by expanding the expressions in powers of r and neglecting 
all but the first terms. Thus it can be shown 

<*»(») 

2 ttT ’ 

^*r) r— *© == H 

(~2)*nla x (x) cos nd 
4vr" +1 ’ 


71=0 

(45a) 

7t>0 




(f) Velocities, vt,. 
Figure 11. — Concluded. 


1 (— 2) n -Tt\a n {x) sin nd 

;. ( P h o)t—*o q 7r ^'^+ 1 ’ 


n>0 


(45b) 


The significance of figure 11 with regard to practical appli- 
cations is more or less obvious. The first step in its use is to 
find the effective length of the cancellation-multipole dis- 
tributions. Since the wing is given, the streamwise varia t ion 
of the cancellation multipoles can be calculated. Actually 
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this variation will extend between the apexes of the enclosing 
Mach forecone and aftercone, a distance of L 0 -\-L' 0 (see 
fig. 10). However, depending on the wing plan form and 
section, the effective lengths of the distributions (the interval 
of principal variation) can be considerably less as illustrated 
in figure 12. Designate this effective length as 2 L, and the 



Fioure 12. — Effective length of multipole distribution. 


distance to the vicinity of the body surface as r„ and one can 
define the parameter r t thus 

(46) 

Using figure 11 and the parameter r e , one can now estimate 
the error incurred by the use of various approximate methods 
for calculating the body shape. A convenient way to carry 
out these estimations is to study the magnitude of the first 
crest of the waves shown in figures 11 (b) through 11 (f), and 
the distance this crest lies from the foremost Mach cone. 
Graphs showing the variations of these quantities with r are 
given in figures 13, 14, and 15. 

By means of the above concepts, let us study briefly four . 
different approximate methods that can be used to calculate 
a body shape. 

SLBNDER-BODY THEORY 

Slender-body theory is represented in figure 14 by the 
straight lines having the slopes, on the log-log scale, equal 
to — (n + 1) where n is the order of the multipole. Since this 
theory amounts to an expansion of the equations for the 
velocities in powers of ?, it obviously represents a good 
approximation when r, is sufficiently small. Notice that for 
a given percentage error the limiting value of r e for which the 
method applies increases as the order of the multipoles in- 
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X 





(c) Tangential velocities induced by second-order multipoles. 
Figure 13. — Continued. 



(d) Radial velocities induced by fourth-order multipoles. 
Fionas 13. — Continued. 


creases. For example, when r,=0.2, <p„ r as given by slender- 
body theory is 19 percent less than that given by exact 
linearized theory for the case shown, whereas <p? r is only 3 
percent less. Correspondingly, the positions of the wave 
crests follow the path predicted by slender-body theory to 
larger values of r as the order of the multipoles increases. 
The latter trend is illustrated by figure 15. 


If for a particular problem r, is small enough for slender- 
body theory to be considered a good approximation, the 
equation for the body shape, r=R{x,0), corresponding to 
the combined wing and optimum cancellation multipoles 
defined in equation (36) is determined by the expressions 
(using equations (45), (36), and (14) together with equa- 
tion (42)) 


dO r <po{x,R,e)' l “ 

dx L U 0 R? _L„ o 


dR ffv(r, Rfi") j 

dx _ Uq — wing 0 


— ) sin 7 id 


4tt 2 R n+1 


2\ n 

- COS 710 


4idR n+1 


r*r ' 

Jo J -L 0 
Jo J 


dx x {x—x^ n cos mpS a "(xi,f) 


/ dx 1 (x—x 1 ) n cos npS w "(?\$) 


(47a) 


(47b) 




(a) <Po r 

Figure 14, — Attenuation of first crest of waves created by multipole 
distributions shown in figure 11. 


(o) <p* r 

Figure 14. — Continued. 
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(d) 

Figure 14. — Continued. 


Approximate methods for finding (< p«IU 0 R) aiK g and (<pt/U 0 ) wing, 
the velocities induced by the wing, can often be used also; 
but these apply to individual cases and cannot be discussed 
here. 

THEORY IFOR LARGE r e 

The asymptotic values for magnitude and position of the 
first wave crest obtained by placing equations (44) into 



(e) <pi, 

Figure 14. — Concluded. 


equations (21) are also shown in figures 13, 14, and 15. 
For to <4 it is clear that this theory can be used when r, 
is greater than about 2. 

CONTROL-SURFACE THEORY 

The appro xim ations inherent in ordinary control-surface 
theory can also be estimated by inspecting figures 14 and 
15, where by control-surface theory one means that the 




THE REDUCTION OE "WAVE DRAG OP WING AND BODY COMBINATIONS 


725 





/ 

/ 

s 

/ 

1 1 

— Slender-body the 

1 1 

K)ry 





/ 

/ 

r—r 1 


r Lar 

/ 

ge 7 the 

ory 





/ 

✓ 

✓ 



















■ 

■ 

■ 

■ 


■ 

■ 

■ 

■ 

■ 

■ 




■ 

■ 


■ 


■ 

■ 

■ 

■ 


■ 

■ 

■ 

■ 

■ 

■ 

Ca) 



■ 



■ 



■ 


■2k 1 1 1 1 1 1 I I I 1 

0 .2 .4 .6 .8 LO 


r 

(a) Radial velocities induced by sources. 

Figure 16. — Position of first crest of waves created by multipole 
distributions shown in figure 11. 



(b) Radial velocities induced by second-order multipoles. 

(c) Tangential velocities induced by second-order multipoles. 

Figure 36. — Continued. 



(d) Radial velocities induced by fourth-order multipoles. 

(e) Tangential velocities induced by fourth-order multipoles. 

Figure 16. — Concluded. 


exact linearized theory is used to evaluate induced velocities 
along a given surface and these values are assumed constant 
for all r in the vicinity of the surface. As shown in 
figure 16, this amounts to assuming <p, and %=» are given 
by a straight horizontal line in figure 14 and by straight 
lines with a unit negative slope in figure 15. Obviously, 
the error in the body shape calculated by this theory increases 
as the amplitudes of the disturbing multipoles increase and 
as the radius of the control surface diminishes. 

One of the simplest applications of control-surface theory 
arises in the study of quasi-cyhndrical bodies. In such 
cases the expression for the body surface can he derived 
immediately from equation (42). Thus, if the amplitudes 
of the cancellation multipoles are small enough and R c , 
the radius of the control surface, is large enough for control- 
surface theory to be considered a good approximation, the 
body shape, r—R(x,6), corresponding to the combined wing 
and optimum cancellation multipoles is determined by 
using equations (19c), (16b), and (36) together with equa- 
tion (42) 


dR r <Pr(x,R e ,d) -\ , ^ <r,CQS7 id p W* $ cosh ( ncosh 1 PR P' 
dx l U 0 4tt 3 R e J_ v VCs-^-W Jo 


$»"(£, ^) cos nip dip 


(48) 


where <t,= 1 for n = 0 and o-*=2 for a)>0. 
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Figure 16. — Comparison of velocity fields given by exaet-linearized 
and control-surface theory. 


A study of optimum fuselage shapes using control-surface 
theory has been carried out by Nielsen (ref. 10) for a constant- 
chor d swep tback wing having a biconvex section and a sonic 
leading edge. The set of interfering singularities used in 
reference 10 are equivalent (the singularities are limited to 
the x axis) to the multipoles used herein. The fuselage 
shapes calculated by Nielsen are thus the same — within the 
accuracy of control-surface theory — as those given by equa- 
tion (48). 


MODIFIED CONTROL-SURFACE THEORY 

A method of modifying control-surface theory to increase 
its accuracy is illustrated in figure 17. Induced velocities 
computed by this method are based on those calculated 
along a given control surface but are extended away from 
this surface by varying their magnitude as (r) r » where the 
value of y n is fixed by the slope of the curves in figure 14 at 
r=r«, r, being defined by equation (46) (fig. 17). With this 
modification equation (48) becomes 

dR r <Pr(?,Rc,0)~\ , ^ ( Rc Y* ** COS 716 

dx~ L U a .W'VW 4 *>R C 


r t-PR, 

J-w 


{x— Q cosh (n cosh 1 d£ 


J SJ'itt'p) cos n\p (49) 

which can be solved using numerical techniques. 

A further refinement of equation (49) can be obtained if 
the position of the induced velocities is also varied according 
to the slope (again at r=r,) of the curves in figure 15. De- 
fining this slope as 5„ see figure 17, and x s as 

x s =x—p(l+8 H )(R—R c ) 


one can see this refinement simply amounts to replacing 
the value of x in the right-hand term of equation (49) by x s . 

ILLUSTRATIVE EXAMPLE— ELLIPTIC WING 

In order that one may be able to assess the practical sig- 
nificance of the preceding sections, the concepts presented 
therein will now be applied to the solution of a particular 
problem. For the basic wing plan form in this particular 



Figure 17. — Comparison of velocity fields given by exact-llnearized 
and modified control-surface theory. 
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example an ellipse will be chosen. There are two good 
reasons for this choice; first, the ellipse is effectively un- 
swopt and places a severe test on the role of body interfer- 
ence in reducing the wave drag at a supersonic speed, and, 
second, for a given volume, the optimum section (i. e., the 
one yielding minimum wave drag) for these wings when 
considered separately has been discovered (see ref. 11) so 
the reduction in wave drag brought about by the body 
will reduce the minimum value possible for such wings when 
flying alone. The drag reductions for the first few cancella- 
tion-multipole distributions will be calculated and compared 
with the total drag of the wing alone, the wing mounted on 
an infinite cylinder, and the wing mounted on a basic body 
of revolution. Finally the details of calculating a body 
shape simulated by the wing source sheet, a source line 
representing a basic body of revolution, and the first two 
optimum cancellation-multipole distributions will be carried 
out. 

THE ELLIPTIC WING 

Consider the elliptic lens specified by the equation 



where the thickness, span, and chord are shown in figure 18. 
The streamwise slope of the upper surface is seen to be 

»)--? <«> 

and the total wing plan-form area S and volume V are, 
respectively, ' 


S=mb "I 



Wing plan form 



Root section 

Figure 18. — Definition of parameters used to study elliptic wing. 
43587(5— 57 17 



Figure 19. — Elliptic wing in coordinates. 


Wave drag. — The wave drag of the elliptic wing repre- 
sented by equation (50) can be calculated by means of 
equation (30) in which, since one wishes now to find the 
wave drag of the wing alone, the a n ’s are set equal to zero. 
The value of S w '(x,6 ) follows by placing equation (51) into 
equation (28) and integrating. Thus 

UoSJ(x,d)=2U 0 cos n ‘ (£+7) sin n)dv= 

2tU 0 cos ii ( . -i? sin p,Vi 

— (^+“2 ~k 

where, by referring the equation of the plan form to the 
coordinates (see eq. (26)) and solving for the points where 
the straight line f= constant intersects the wing edges, one 
finds — see figure 19 


6 2 gsin ^dcab-yjg} cos 2 sin 2 /i— £ 2 cos 2 p 

M ar cos* sin 2 p 


Hence, 


X / (^,e)=- (a , + ff a c 5 083 ey ^+bW COS 2 g-2 2 (53) 


From the relation 


Z 2 (0)=a 2 +6 2 |3 2 cos 2 0 (54) 

the wave drag can be expressed in the form (integrating once 
by parts) 


d_ i r*- M C L(e > f £ »> r4fan : 

2 47T 2 Jo J-L(0) ( ^ l J-L(e) ^ l_h*(0)J 

YL\0)-2tn 


L\e)~ 2&n &y£*(g)-k* 

JWe Fi?J fi-fe 


Further integration yields 


Jo (a 2 + 6 2 j9 J cos 2 BY 
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Values of <£y/££ 

Point Exact Ackeret 
number linearized wave — 
theory t{x-r)/a z 

1 .044 .050 

2 -.053 -.050 

3 -.102 -.100 

4 -.104 -.100 



Figubb 20. — Points at which exact-linearized theory was compared 
with an Ackeret wave in plane of symmetry for elliptic wing. 


Finally, tlie wave drag can be expressed in coefficient form, 
based on the total wing area irab, as 



Equation (55) represents the lowest value of wave drag 
possible for a wing having an elbptic plan form and fixed 
volume. This equation was first derived by Jones in 
reference 11. 

The velocities induced by the wing source sheet in the 
vicinity of the fuselage. — Later, when one wishes to calcu- 
late a stream surface in the presence of the source sheet that 
simulates the wing given by equation (50), it is necessary to 
know the velocities induced by these sources at the body 
surface. Hence, the value of <p T induced by the source sheet 
was calculated at the four points indicated in figure 20. 
As it turns out, these values are so close (see the figure for a 
numerical comparison) to those obtained by assuming the 
source sheet to he two-dimensional with a chordwise intensity 
identical to that along the root section of the elliptical sheet 
(i. e., using the Ackeret wave generated by the root section) 
that the effect of the wing can be assumed to be given every- 
where in the vicinity of the body by the latter velocity field 
if (as will be the case in subsequent application) the surface 
of the body passes through the region shaded in the figure. 
That is, the effect of the wing in the equations for the fuselage 
shape (such as eqs. (47), (48), or (49)) is assumed to be 


1 t , . . . . — forO<0<w 

U*= T a* {xTpr sm*)sm0, +for7r < 0 < 2 ^ 

1 f , „ . „ — forO<(0-<7r 

“n« )co8 «, + f O rw<0<2w 


THE OPTIMUM CANCELLATION MULTEPOLES 

One can now find the strengths of the multipoles along the 
x axis which induce around a cylinder of infini te radius a 
momentum field identical to that created there by the elliptic 
wing. The negatives of these variations are, according to 


equation (36), the optimum cancellation multipoles. Hence, 
combining equations (53) and (36) 


(%) = 


2tabUoX f 2T V® s + cos 3 8— x 1 * 
* n <-j9) n Jo (o. J +& J 0 s cos 'ey 


cos nO do 


(67) 


where <r n =l for n= 0 and <r»=2 for n> 0. Particular vari- 
ations of a%°(z) are shown in figure 21. These results aro 



(a) n=0 



Figure 21. — Variation of nth derivative of nth-order cancellation 
multipoles for elliptic wing. 
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for «=0, 2, and 4, since a™ (x) for any odd n is zero by 
symmetry, and apply when the wing plan for and free- 
stream Mach number are related by 


a 4 
b/3~3^ 


(58) 


which contains the particular- case for which the Mach 
number is ^2 and the aspect ratio is 3. It is apparent that 
there are at least -n+1 roots to a^(x) for —L 0 <z<L 0 . 
This follows immediately from equation (38) and is true in 
general. As a result the curves for the higher values of n 



X 


(a) »=2 



(b) n=4 

Figtjbis 22. — Variation of nth-order cancellation multipoles for elliptic 

wing. 


become increasingly wavy and, correspondingly, increas- 
ingly difficult to evaluate numerically. 

Figure 22 presents the values of a n (x) for the same elliptic- 
wing Mach number relation given by equation (58). Notice 
that each of these curves has only one root (they necessarily 
have at least one) in the interval — L„<^x<C.L 0 and is in- 
creasingly smooth with increasing n. The latter follows 
from equation (37) and the fact that the first, n derivatives 
of these curves must, in general, be continuous. For exam- 
ple at x= ±L 0 the first four derivatives of a^x) must vanish. 

Wave drag. — One can now calculate how much the wing- 
alone drag is reduced when combined with each successive 
opt imum cancellation-multipole distribution. If D % denotes 
the drag saved by the nth-order cancellation multipoles, 
then by equation (35) 

~=~l TU?j- L dXl ^- L dx ' 2a ^ + ' ) + ” (■ r -) 1 (59) 

where L 0 is the maximum value of L(O) as given by equation 
(54) 

L^at+b*? (60) 

The total drag saved by means of the first m multipolo dis- 
tributions, would, by equation (34), be 


m 

D=2JJ 0 -\-'^2iD n (61) 

Using equations (53) and (36) to define the ai n+1) (r) in 
equation(59), reversing the order of integration, and inte- 
grating once by parts, one finds 
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It is apparent from figure 23 that this can bo written 


D n 16 (Jab)’ 


or 
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The total drag saved by using all the cancellation multipoles 
is, by definition, 
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Figure 23. — Regions used in developing equation (62). 


and since 

x=t— S - sin 2nx cos 2nx 
4 n 

this is equal to the drag of the wing alone, as it, of course, 
should be. 

The reduction in wave drag as the wing is combined with 
the first three optimum cancellation multipoles is presented 
in figure 24. In studying figure 24, one sees, as the Mach 
number approaches 1 (i. e., /3 — >0), more and more of the 
original wing wave drag is destroyed by a line of simple 
sources alone. Further, the value of r, which can be writ- 
ten — see equations (46), (54), and figures 21 and 22 — 


* Va s +6 2 0 2 


(63) 


tends (for a fixed average distance to the body surface r e ) 
to zero as the Mach number approaches 1 ; and this, in turn, 
means that as /9 goes to zero the effect of the multipole 
strengths on the body shape can he calculated using slender- 
body theory. 



Figure 24. — Portions of elliptic-wing wave drag oreated by various 
equivalent multipole distributions. 


. 10 , 



x 

(a) Radial velocities induced by sources. 

Figure 25. — Velocities induced by the elliptic-wing cancellation 
multipoles at the control surface where pl?/L.= 0.148. 

When various orders of multipoles are distributed along 
a line, one can show the cross-sectional area normal to the 
free stream of the simulated body as given by slendor-body 
theory is a function of the source distribution only (soe 
Appendix B). Coupled with the discussion in the preceding 
paragraph, this can be used to demonstrate that, for Mach 
numbers close to 1, the “supersonic area rule” proposed by 
Jones (ref. 12) and Whitcomb and discussed in reference 13 
gives a good approximation for the wave drag of an elliptic 
wing and body combination which is symmetrical with 
respect to the plane of the wing. 

The induced velocity field. — A method for calculating the 
velocity field induced by the multipoles when a n (x) is given 
numerically is presented in Appendix C. By means of this 
method, velocities induced by the ao and a 2 multipole dis- 
tributions shown in figure 22 have been calculated for r 
equal to 0.148 and the results are shown in figure 25. 
Since the distributions in figure 22 were for the particular 
case a/6i3=4/3ir, it is evident from equation (63) that the 
values in figure 25 apply to the case r,/6 equal to 0.161; 
that is, when the body radius is about 16 percent of the 
wing semispan. 

For comparative purposes, the values given by slender- 
body theory are also shown in figure 25. The degree of 
agreement between the two curves is consistent with the 
results shown in figures 14 and 15. 

INTERPRETATION OF DRAG REDUCTIONS 

Comparison with wing mounted on a circular oylinder. — 
With regard to figure 24 one should be careful to notice 
that the drag of the wing alone has been used for the 
reference drag. The drag reductions shown, therefore, repre- 
sent gains brought about by interfering with the velocity 
field induced by a planar source sheet, or, in terms of a 
combination with an upstream cylindrical stream surface, 
gains made by modifying a body, shown in figure 26, which 
bulges behind the wing leading-edge Mach wave in accord- 
ance with the velocities induced there by the source sheet. 
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(b) Radial velocities induced by second-order multipoles. 
Figure 25. — Continued. 


Obviously, from this viewpoint, a considerable reduction in 
drag con be brought about merely by eliminating the bulge, 
thereby making the body a circular cylinder throughout. 
Mathematically, such a procedure amounts to using a certain 
set 0 of cancellation multipoles along the x axis behind the 
point — L„, and, if the drag of this resulting combination 
were used as a reference, the gains shown in figure 24 would 
be diminished. 

An approximate way to estimate the drag of a wing 
mounted on a circular cylinder is illustrated in figure 27 
and consists merely of subtracting from the wing source 
sheet those sources blanketed by the body. Using the 
subscripts 1, 2, and 3 to designate the wave drags of the 
individual wings as indicated in figure 27, Jones (ref. 15) 
has shown that if wing 1 is an elliptic wing with a bi- 
convex section and wing 2 lies entirely within the plan 
form of wing 1, then 

U 3 -U 1 (l+^)+U s (64) 



Fiourh 20. — Wing and stream tube simulated by planar sheet of sources. 

• The exact evaluation of mnlHpole distributions necessary to simulate a circular cylinder 
for the entire body length has been studied In reference 14. 



X 

(c) Tangential velocities induced by second-order multipoles. 
Figure 25. — Concluded. 


where V 2 and Vi are the volumes of wings 1 and 2, respec- 
tively. 

For a supersonic Mach number, D 2 is closely approximated 
by the wave drag of a rectangular wing having the same 
section and aspect ratio. If A 3 , 4 Ra, and r t are, respectively, 
the aspect ratio, plan-form area, and thickness ratio of the 
rectangular wing, its drag can be expressed in the form 





Equivalent wing-body combinations 


Figure 27. — Method of approximating wing mounted on a cylinder 
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Further, if is the aspect ratio of the elliptic wing, one can 
show — see figure 27 and equation (52) — 

(67) 


The drag of the elliptic "wing follows from equation (55) and 
can be ■written 

^=4 abN t (68) 

To 

where 


Ni=v 


1+2 

J 4 N 

WAi/ 

•» 

>' 

X 

4 Vl v - 

*Mi/ J 


(69) 


Finally, therefore, equation (64) can be put in the form 



and the ratio N 2 /Ni is a function of the parameters R/b and 
pA l only. 

By means of equation (70), the dashed curves shown in 
figure 28— 



Figure 28. — Drag of various multipole distributions compared with 
drag of wing mounted on circular cylinder. 

representing approximately the wave drag of a wing mounted 
centrally on a circular cylinder — were calculated. Though 
considerable drag reduction is indicated by adding just 
those multipoles necessary to make the body cylindrical, 
it is apparent the total wave 'drag can be reduced further, 
for the range of parameters shown, by using only the first 
two optimum cancellation-multipole distributions, a<,(a;) 
and aj(r), given by equation (36). 

Comparison with wing mounted on a basic body of revo- 
lution. — Figure 28 shows the effect on the wave drag of 
adding the op timum cancellation multipoles either to the 
wing alone or to the combination of an infinite circular 
cylinder and a centrally mounted wing. Estimates of their 
effect when added to a wing mounted on a basic body of 
revolution can also be carried out. In order to present these 


estimates, however, the results of the following two theorems 
due to R. T. Jones (ref. 15) are needed. 

1. Designate the closed body of revolution which, by 
slender-body theory, has a minimum drag for a fixed volume 
and length as a Sears-Haack body. Then the total wave 
drag of a Sears-Haack body and any other body of revolution 
or any centrally mounted thin wing which he entirely within 
the Sears-Haack body’s enclosing Mach forecone and after- 
cone is given by the equation 

D=D SB (l+y^)+D 2 (71) 

where: 

D sa wave drag of Sears-Haack body alone 
D 2 wave drag of other body or (exposed) wing alone 
Vss volume of Sears-Haack body 
V 2 volume of other body or (exposed) wing 

2. Designate the body of revolution which, by slender-body 
theory, has a minimum drag for a fixed base diameter and 
length as a Kkrm&n ogive. Then the total wave drag of a 
Ktrmkn ogive and any other slender body of revolution or 
any centrally mounted thin wing which he entirely within 
the ogive’s enclosing Mach forecone and aftercone is given 
by the equation 

D=D X +D a (72) 

where: 

D k wave drag of K&rm&n ogive alone 

wave drag of other body or (exposed) wing alone 

In order that the theoretical results could bo tested by 
wind-tunnel experiments, a basic body of revolution having 
a finite base area was chosen. Such a body can be simulated 
by a combination of the source distributions which produce, 10 
separately, the Sears-Haack body and the Kfirmfm ogive. 
Thus, if 2 1 is to be the body length, the line of sources 

(73) 

simulates (by slender-body theory) a body of revolution 
having a total volume V equal to V aB -\-V K , a cross-sectional 
area given by 

Sl*)-^! [rVF^-W 3 (f+sin- 1 £)]+ 

J^FG’-Z *) 3 ' *, -1<X<1 (74) 

and a base area S(l ) equal to V E /l. 

The wave drag of a wing mounted on this basic, unmodified 
body will now be calculated. Just as was the case in study- 
ing the wing attached to an infinite cylinder, the assumption 
is made that the wave drag of this combination is the same 
as the wave drag on the configuration simulated by super- 
imposing the singularity distributions which create sepn- 

if The source distributions simulate the Sears-Haack body and the KArmAn oglvo only 
when slender-body theory Is used to calculate the body shapes. If linearized theory Is used 
instead, the body shape will, of course, be somewhat different. However, the subsequent 
results and conclusions are by no means limited by the assumptions pertaining to slender- 
body theory. The latter theory is used only to obtain an estimate of the body volumes or 
to study cases for which It gives results that are not significantly different from those given 
by linearized theory. 
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rately the exposed wing panels and the body of revolution. 
With this assumption, the wave drag can be written ex- 
plicitly in terms of the wing and body geometry by applying 
equations (71) and (72). Hence, 

D^Dx+Dsafl+^+Ds (75) 

where D 3 is the drag of the exposed wing panels alone, given 
by equation (70) and shown for various values of R/b in 
figure 28, and V 3 is their volume (see fig. 27) 

f— I t) <76) 

Since D K and D S h, the wave drags of a Kdrinttn ogive and a 
Sears-Haack body flying alone, are well known to be 



the wave drag coefficient of the unmodified combination, 
based on the complete wing area nub, can be expressed as 

D „ TV+8HV , \UV SI ,(n 8 R'\ , 
qmb~^ D - 7 T*l*ab + **1* \2 3 b )' 



where Ni and N 2 are defined in terms of Mach number and 
wing-body geometry in equations (66) and (69). An ex- 
ample of the variation of C D with Mach number for the par- 
ticular combination shown in figure 29 ( R/b was set equal to 
0.181) is given by the dashed line in figure 30. 

It is now possible to find how much the drag of this un- 
modified combination can be reduced by means of the opti- 
mum cancellation-multipole distributions used to derive the 
results shown in figure 24. Again applying equations (71) 
and (72), one can show 

D=DK+D SH (l+^)+D i (79) 



Maximum thickness of wing 
along center line, / =0.234 

Total volume of body = 44.60 
VSH-- 3172 
V K - 12.88 



Figure) 30. — Drag variation for modified and unmodified wing-body 

combinations. . 


where V* is the total volume and D t the total wave drag per- 
taining to the wing (now the complete wing including the 
portion blanketed by the body) and the multipoles. How- 
ever, within the accuracy of the approximation — being, in 
fact, exact within slender-body theory, see Appendix B — 
the volume added by the wing is subtracted from the basic 
body by the optimum cancellation-source distribution so 
that Vi is zero. Further, if Ni is the value of DjD w read 
from figure 24 for a specific value of fiAi and a specific num- 
ber of multipole types, one can readily show 

DfJ^N.Ni (80) 

where N x is defined in equation (69). Hence, the drag of 
the unmodified combination can be reduced to either 


y ' ,+s [>»+ teS (!-| m ■ 




if the same total volume is maintained (maintained, as is 
obvious from an inspection of the equation, by increasing 
the value of the Sears-Haack portion of the basic body an 
amount equal to the volume of the exposed wing) or to 



Fiqure 29. — Dimensions of wing-body combination analyzed. 


if the volume of the fuselage is reduced by an amount equal 
to the wing volume. 
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Actual 

cross section 
at station 



B.5 

I 



I 1 




sides by the cancellation multipoles (see fig. 31), decreasing 
the average body radius in the wing region from about 1.00 
to about 0.89. One must be careful to notice that tho solid 
curves represent minimum (relative to the special method 
being discussed) values which can be obtained by a specific 
design at a specific Mach number and do not represent the 
variation of wave drag with Mach number for any given 
combination. 


THE BODY SHAPE— BTBST CALCULATION 

The final step in studying the effect of the optimum 
cancellation multipoles, defined in equation (36), is to find 
the distorted body shape which they produce in combination 
with the wing and a basic body. The decision was made to 
calculate a body shape which would be optimum at a Mach 
number equal to -J2. The details of tho wing and body 
geometry are given in figure 29 and the basic body parameters 
V K and V s were interpreted in terms of source strength by 
equation (73). 

It was apparent from the results of figure 26 that, for the 
values of r, and b given by figure 29, the velocity field induced 
by the first two optimum multipole distributions can bo 
calculated with good accuracy using slender-body theory. 
Combining the values of <p T and <po so calculated with those 
induced by the wing, given by equation (56), and thoso 
induced by the basic body, using also slender-body theory to 
interpret equation (73), one can find tho body shapo by 
solving the two simultaneous nonlinear differential equations 
presented as equations (43). These were solved numerically 
by the method outlined in Appendix D and the results were, 
unfortunately, unrealistic. Figure 32 shows an example of a 
streamline close to the 0=0 plane and the crossing of such 
streamlines obviously invalidates the solution. 





Note- vertical scale five times horizontal 


Figure 31. — Body shape having favorable interference wave drag at 
flf=1.41 when combined with elliptic wing. 


The results expressed by equations (80) and (81), when 
applied to the first two optimum cancellation-multipole dis- 
tributions, are shown for the geometrical parameters pre- 
sented in figure 29 by the solid curves in figure 30. The 
value of R/b used for the solid curves was 0.161 instead of the 
0.181 value used to calculate the dashed curve. The smaller 
value was used since the modified body is drawn in along the 


BODY SHAPE— SECOND CALCULATION 

The failure observed in the first calculation has a simple 
enough interpretation. For the chosen wing the basic body 
was too small in diameter at the -wing-body juncture to per- 
mit the use of the first two cancellation multipoles in their 
entirety. 

Several avenues of approach are yet available. One could, 
for example, maintain the same -wing and basic body but 
reduce the Mach number, one could start with a larger basic 
body, or one could lower the thickness ratio or aspect ratio 
of the wing, thus diminishing tho strength of the cancellation 
multipoles. All of these, however, are modifications of the 
basic conditions or basic restraints and as soon as such 


/ 



Figure 32. — Streamline when starting body radius is too small. 
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restrictions are abandoned it must be remarked that no 
matter how low the wave drag of a set of nonlifting, volume- 
enclosing surfaces has been, made, another arrangement of 
the same volume within a finite space will give a lower value, 
unless, of course, the wave drag of the first arrangement is 
already zero. Therefore, instead of modifying any of the 
initial restraints, consider the following alternative: 

How much can the drag be reduced by using only a 
portion of the first two cancellation multipoles so that a 
real body would still be simulated? 

In order to answer the above question, examine briefly the 
first calculation. Notice, from figure 25, that the cause of 
the body collapse is attributable to the large values of 
and ipe induced by the second-order cancellation multipoles, 
a 3 (x). Hence, let a 0 (x) be maintained at its full value and 
reductions permitted only in the magnitude of a 3 (x). To 
carry out such a procedure efficiently, one must be able to 
determine the effect of a given variation bf a 2 (x) on the body 
shape. Fortunately, Graham, in reference 16, has developed 
a method by means of which the relation between a 2 (x) and 
body shape can be quickly estimated. Graham has shown 
if 

1. a rectilinear distribution of second-order multipoles 

of strength ce 3 (x) is placed along the x axis (see fig. 

33) in a supersonic stream 

2. slender-body theory is used to evaluate <p T and tpe, 

3. R=R e is the radius of a circularly cylindrical tube 

for — £ 0 >r>— 0 ° , 

4 - S-L.^ )dlh 

then 

^T l+2y*(r) cos 20+ Vl +4F*(3Q cos 20+4[-F*(a;)] a \ 1/ ' > 

(83) 

is the continuation of the stream tube for L 0 . 

Since the initial strengths of the cancellation multipoles 
are negative, F*(x) is negative and the critical value of R 



Figure 33. — Definition of parameters used to study effeot of seooDd- 
order multipoles on fuselage shape. 

430875—07 48 



Figure 34. — Variation of R/R e with — F*„„. 

occurs along the plane 0=0 or tt. The variation of R/R c 
with F* max is given in figure 34 from which one can see that 
the maximum value of |F*| must be less than 0.5 if the 
simulated stream tube behind the plane x——L 0 is to be real. 

The problem can now be continued, using Graham’s 
result as a guide, by assuming- the critical body radius in the 
more complicated source and multipole arrangement is 
principally determined by F*{ 0) — the parameter governing 
the body indentation at the center of the cutout and at the 
wing-body juncture. In the first place, since there is no 
interference between different orders of multipoles, it is 
necessary to consider only the drag produced by the second- 
order multipole. Appendix E presents a method for finding 
the optimum distribution of the second-order cancellation 
multipoles for a given wing and a fixed value of F*( 0). 
The resulting wave drag is given in equation (Ell). At a 
Mach number equal to -y/2 and for the basic wing and body 
parameters presented in figure 29, F w *( 0) (defined by eq. 
(E8)) equals —2.90 and the reduction in DilgS, tho amount 
of drag caused by the wing second-order multipoles alone — 
see equation (62) and figure 24 — is shown in figure 35 for a 
range of F*( 0). Variations on the strength of various com- 
binations of second-order multipoles are shown in figure 36. 
It is important to notice that for a given percentage reduction 
in the maximum strength of the multipoles the resulting 
percentage reduction in D 2 is much larger. 

The strengths of a 2 (x) shown in figure 36 must now be 
combined with the zero-order multipoles and wing source 
sheet, and the combined velocity field used to calculate the 
shape of the new body. Using again slender-body theory to 
evaluate the velocity field induced by the multipoles and the 
numerical methods given in Appendix D to compute tho 
streamlines, one finds, by restricting the distribution and 
strength of the second-order cancellation multipoles to their 
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Figure 35. — Reduction in drag caused by second-order multipoles. 



Figure 36. — Variation of second-order multipoles for several values 

of F*(0). 


optimum values corresponding to the restraint 11 F*( 0) = 
—0.6, that a real as well as reasonable body shape results. 
The details of this shape are presented in figure 31 and then- 
general interpretation is discussed in the next section. 
Finally, using the value ^(0) = —0.6, the drag curves shown 
in figure 30 were reinterpreted, and the results— which rep- 

u Figure 34 gives 0.5 as the marimnm permissible value of — F*(x). However, that value 
Is based on a distribution of ca(?) alone, and In our more general case the added velocity 
field caused by the presence of the other singularities permits the larger value. 



Figure 37. — Drag of unmodified and modified combinations with real 

body shapes. 

resent an estimate of the amount the wave drag of an 
elliptic wing mounted on a basic body of revolution can be 
reduced by realistic body distortions — are shown in figure 37. 

discussion of results 

It is possible to gain some insight into the reasons for the 
various body distortions shown in figure 31 by inspecting, 
in another light, the body shape first calculated. Consider 
first the elliptic wing at the top of figure 38. The air over 
the forward part of this wing, when it is alone in a super- 
sonic stream, is compressed (mathematically, the sign of <p z 
is negative), the compression being greatest near the leading 
edge. On the other hand, the air over the after portion of 
the wing is undergoing an expansion, the magnitude of which 
is greatest near the trailing edge. Consider now, in combi- 
nation with this wing, a body which is to have a shape 
providing favorable interference. It is apparent that the 
body should cast expansion waves over the forward portion 
of the wing, destroying the compression there, and absorb 
the expansion waves coming from the wing afterportion. 
Or in another fight, the positive pressure on the forward 
region of the wing (one can use the equation C p ={p—p 0 )l<L — 
—2<p x /U 0 for the pressure coefficient) should be reduced as far 
as possible by a wave shed from the body and having large 
negative pressures where it comes in contact with the wing 
forward region. 

Since waves in a supersonic flow field are fundamentally 
associated with the slope of the disturbing surface, the afore- 
mentioned favorable interference fields would be created by 
a body having, longitu dinall y along its surface, slopes such 
as those shown in the lower part of figure 38. This is ex- 
actly what the solution obtained from the calculation of the 
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Figure 38. — Regions of expansion and compression on wing with 
unmodified and modified fuselage. 

first body shape tried to establish since the fuselage near the 
plane of the wing (the portion most strongly affecting and 
being affected by the pressures on the wing) and ahead of the 
wing chordwise center line was distorted in a manner that 
caused an expansion across the wing entire forward portion. 
The difficulty arose because the fuselage was not wide enough 


to provide the longitudinal extent of favorable slopes neces- 
sary to create the positive pressure called for bj r the wing 
forward compression region, and the body streamline near 
the wing root, following a path such as that shown by the 
line in figure 38, crossed the body center line bofore it reached 
the wing chordwise center. 

Consider now the second body calculated in the previous 
section. In this case an additional restraint was imposed 
which, effectively, 'fixed the maximum body indentation. 
Subject to such a condition, an optimum interference field 
was discovered. If the resulting fuselage shape is inspected 
near the plane of the wing, surface slopes are found similar 
to those shown in figure 39. The following discussion is 
intended to show that, from a physical viewpoint, this ar- 
rangement is reasonable. 

Most of the wing pressure drag occurs on the -wing inboard 
portions. Hence, for a fixed maximum fuselage indentation, 
it is beneficial, from an over-all point of view, to create 
initially a compression wave, which increases the pressure 
drag on the forward portion of the wing tip but provides a 



Figure 39. — Regions of expansion and compression on wing with 
fuselage distortion shown. 
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succeeding extent of fuselage having slopes that generate a 
strong expansion wave over the forward portion of the wing 
inboard section. Similarly, the final portion of the body is 
forced to have a region of unfavorable interferences where 
the expansion waves from the wing tips combine with body 
expansion waves to increase the local drag (i. e., increase the 
local suction pressure) in order that the over-all interference 
effects are as beneficial, under the given restraints, as possible. 
This arrangement (i. e., unfavorable interference near the 
wing tip and favorable interference near the wing root) is 
given further support by the attenuation property inherent 
in three-dimensional waves. Thus the pressures induced 
by the body on the wing tips are not as strong, for a given 
generating surface slope, as those induced on the inner por- 
tion of the wing, simply because the tips are farther from the 
disturbing surface. 

Although these considerations are somewhat oversimplified 
(the shape of the upper part of the body has been completely 
ignored in estimating the effect of the waves), the longitudinal 
variation of surface slopes near the plane of the wing and the 
resulting body streamlines there are, from a physical point 
of view, reasonable. 

In order to support the above conclusions, the source and 
multipole distributions simulating the final modified body 
shown in figure 31 were used to calculate (see Appendix C) 
u/U a in the plane of the wing near the root section. The 
values of uJU 0 induced by the wing sources along these 
sections were assumed to be the same as those induced by a 
two-dimensional biconvex section having the same local 
chord; that is, tip effects were neglected. These values for 
body and wing were added and the resulting pressure distri- 
bution, shown in figure 40 {O v =— 2 u/U 0 ), were obtained. 



Figure 40. — Variation of pressure and section-drag coefficient on 
sections indicated. 


The results are similar to the estimates presented in figuro 
39. The large drag saving near the root section is illustrated 
in figure 40 by the graph showing the low values of section 
drag coefficient along the inner portion of the tying. 

Another important characteristic of wing-body combina- 
tions designed to have low wave drag is also illustrated in 
figure 40. As shown in the graph of u/U 0 , over the surface 
of a two-dimensional biconvex section the air is everywhere 
accelerating in the streamwise direction. In studies con- 
cerning the effects of viscosity on the fluid flow and, in 
particular, studies concerning the boundary layer, this posi- 
tive fluid acceleration is referred to as a favorable pressure 
gradient. If the flow is laminar in the vicinity of the leading 
edge of a smooth wing and the pressure gradient is every- 
where favorable, the flow tends to remain laminar and 
unseparated over most of the wing chord. Notice that the 
modified wing-body combination has a line of zero pressure 
gradient extending along a Mach line downstream from a 
point near the body and wing leading-edge juncture. Im- 
mediately behind this line the pressure gradient is unfavor- 
able which gives rise to the possibility of flow separation or, 
at least, transition from laminar to turbulent flow there. 

COMPARISON WITH EXPERIMENT 

The modified wing-body combination shown in figuro 31 
was tested in the Ames 2- by 2-foot transonic wind tunnel. 
The Reynolds number of the test, based on the mean aero- 
dynamic chord, was approximately 1.5 X10 6 . This com- 
bination had an exposed wing volume of 3.44 cubic inches 
and a body volume equal to 44.60 cubic inches, for a total 
volume of 48.04 cubic inches. As a control, an unmodified 
combination composed of the same elliptic wing mounted 
on a body of revolution (the area distribution of which was 
determined from equation (74) with /=10.5, V K = 12.88, 
and V sb =29.02 cubic inches) was tested. The exposed 
wing area in the unmodified combination was 3.32 cubic 
inches and the body volume was 41.00 cubic inches, for a 
total volume of 45.22 cubic inches. Thus, the unmodified 
combination had the same body length as the modified one 
but less volume. 

The wave drag at M=1.41 of the combination shown in 
figure 31 has already been calculated and presented hi figuro 
37 by the curve pertaining to real body shapes. By use, in 
equation (78), of the values of V K and Vs,, mentioned above 
and a value of 0.176 for R/b, the wave drag for the unmodified 
body was calculated throughout a supersonic Mach number 
range. The theoretical results obtained, for body configura- 
tions are shown by the dashed curves in figuro 41. 

The wind-tunnel results for the total drag on both con- 
figurations are shown in figure 42 for 0.6<M<1.4. Notice 
that three groups of data are shown. The lower one repre- 
sents the unmodified body alone, the middle one represents 
the modified and unmodified combinations with no fixed 
transition, and the upper one represents both combinations 
with transition fixed along the leading edge. The models 
tested with natural transition did not show the predicted 
drag reduction. As was pointed out in the discussion of 
figure 40, however, the adverse pressure gradients on the 
modified model could be inducing transition in the vicinity of 
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Fioube 41. — Variation with Mach number of theoretically and experi- 
mentally determined drag coefficient. 


the 45° line behind the juncture of the wing leading edge with 
tho fuselage. This, in turn, would cause the wing of the 
modified model to have a larger area covered with a tur- 
bulent boundary layer and, hence, cause the drag of the 
model to increase. In order to separate the potential and 
viscous effects, the transition-fixed tests were made. If the 
experimental wave drag is taken to be the difference between 
the drag at a supersonic Mach number and the drag, at 



Figure 42. — Experimental values for drag coefficient. 


M= 0.6, the resulting values of experimental wave drag are 
as shown in figure 41. 

Figure 41 shows that the experimental reduction in 
wave drag brought about by the modification agrees with 
that predicted by theory. Both theory and experiment 
show a reduction of about 0.0015 in the drag coefficient at 
the design Mach number (1.41), and the experiment further 
shows an average reduction of 0.0020 over the Mach number 
range 1.2<M<1.4. A further study of figure 42 shows that 
the difference between the experimental and theoretical 
wave drags shown in figure 41 for the wing-body combina- 
tions is nearly the same as the difference between experiment 
and theory for the body alone. 

Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 16, 1956 




For convenience, take the normalized form of the wave 
equation in Cartesian coordinates, thus 


i) 2 tp d*<p 

d^~dy~bF- 

and define the Laplace transform of <p(x,y,z) by 


v(s,y,z)=J o 


APPENDIX A 

DERIVATION OF THE OPERATIONAL FORM OF THE WAVE EQUATION 

Further, integrating by parts gives 

sf). 

Hence, 

jx3-W)*-*3-s i- 


e a ^(x,y,z)dx 


(Al) 


(A2) 


Now if x=f(y,z ) is the equation of the foremost Mach cone 
or Mach cone envelope and f(y,z) > 0, it is apparent 

r v (tvz)e-*dx= r 

W ^ X ’ y ’ Z)e J , e dy^ Zyy e \ZyyLt 

(A3) 

since is, but (j^J is not, necessarily zero. From 

equation (A3) we see 


'z-r 


(A5) 


(A6) 


\dx~^by Z>y'Z)z Z>z/ X -/ 
(A7) 

The last term on the right is the directional derivative of 
the perturbation potential along the surface x=J(y,z). 
This is, of course, along the so-called conormal. Since <p is 
a constant on the forward envelope, its gradient along the 
envelope is zero and 


APPENDIX B 

ON THE VOLUME OF BODIES CALCULATED USING SLENDER-BODY THEORY 

multiply by dd, and integrate 


The following proof shows that in a rectilinear distribution 
of singularities, only the sources contribute to the total cross- 
sectional area of the simulated body and, hence, to its volume. 

According to slender-body theory, the velocities induced 
in the field by distributions of multipoles along the x axis can 
be written 




ay(x) 

2irr 

(—2Yn\a n (x) cos nd 


4 vr' 


,a+l 


71=0 

7l>0 


1 (—2)*n\a n (x) sin nd 

r v *»~ 47T7--+ 1 


(Bla) 


(Bib) 


Further, we have derived — see equation (41) — neglecting 
only second-order effects, the equation representing the 
boundary condition for the body, thus 


p. ® 


0* T 2 2# 

Combine equations (Bl) and (B2) 

0 TT ~r\ Z)R “ f cos nB wind Z>R\ 

2 * U ° E ( 2) n nla n (x) -^r ^ ) 


(B2) 


2 ^4Jo t^= 2 ^)+ 

S (— , Z) n n\a n {x) f 
i Jo 


R cos nd — sm nd ^ 


~W+ r 


7 


or 


2 *U„ ^ S(x)=2*a 0 (x)+± (~2) n (n—l)\a n (x) J" f Q dd 

(B3) 

Since the integrand in equation (B3) is a periodic function 
in 6, we have 

bS_ a v (x) j-n .. 

bx~ U 0 K ) 

which shows the simulated-body normal cross-sectional area 
to be dependent only on the source strength. Further, the 
total volume is given by 

V= J' S(x)dx=J ! ' {l'-x)S , (x)dx+(l'+l)S(-l) 

and when S(l')=S(—l)=0, there results 

V=— ^-J* xa 0 (x)dx (B5) 
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APPENDIX C 

ON THE CALCULATION OF VELOCITIES INDUCED BY ARBITRARY SOURCE DISTRIBUTIONS 

Then 


The potential and velocity fields represented by equations 
(18) and (19) are difficult to evaluate analytically even if 
xl„(£) is a simple function. However, the calculations can 
be reduced to a relatively simple process. First, let equation 
(18) be expressed in terms of the dimensionless variables 
x,£, and 7 where 

x=x/L, 


7=pr/L 0 


(Cl) 


/=-« i n- ? Aj&te 
v '®’ *J-. jraa 


<p>(x,r,ey 


VOMJ* 

cos 2 e p- 7 A 2 (i)[2(5-|) 2 -p]dF 
2 t r 


if; 


W(x-D 2 -p 


(C2a) 


(C2b) 


and so forth. Consider next the variation of A B (|) shown 
in figure 43 and represented by the equation 


A»(D= 


( (190a-360m)+(3907i-760ra)£+(200a-400mjF, -1<£<— 0.9 

n , 


-0.9<£ 


(03) 


As seen in the figure, A n vanishes at £= — 1, is a parabola 
between —1 and —0.9 (assuming the values m at £= — 0.95 
and 72 at £= — 0.9), and the straight fine, A„(J)=n, for £ 


greater than —0.9. The velocities induced by a multipole 
distribution given by equation (C3) can be calculated in a 
straightforward manner in the two regions — l+7<x< 
— 0.9+7 and — 0.9+7 <x. For example, if 


vrMcffl J 

200 =H( 


-(z+0.9)V(7+1-0) 2 -P+P1 n ^jJ+Vg+l-O) 2 r* f -l+?<x<-0.9+7 


-(x+0.9)V(x+1.0) 2 -P+(x+1.0)V(x+0.9) 2 -P+Pl^ 


x+1.0+V(x+1.0) 2 — r 


(C4) 


X+0.9+ Vtx+0.9)* — r 


, — 0.9+r<x 


irrN^ix) 

100 


(x+0.95)V(x+1.0) 2 — P — Pin. a;+1,0+ Vg+1-0 Y 7 2 j 9+? 


= -< 


(x+0.95)V(x+1.0) 2 -P-(x+1.05)V(7+0.9) 2 -P-Plu ;| ^ j ^ -0.9+7<x 

X+0.9+ V(7+0.9) 2 — P 


(Co) 


then <p rr cap be written 

^07=mM Of (x)+7iA r o? (x) (C6) 

Now, if one is given a distribution of sources that is 
composed of, or is adequately approximated by, a series of 
20 equally spaced parabolic arcs, equation (C6) can be used 
for each individual arc and the results superimposed for the 
complete solution. To this effect, define m t and n { in terms 
of A 0 (x) by 

y (P7) 


so they represent the magnitudes shown in figure 44. Then, 
if [x] denotes the greatest interger contained in x(e. g., 
[6.34] equals 6), the equation for the radial velocity becomes 


Vr ~ § [ m ' M ‘ 7 G _ Tr) +7ltAr ^( 2_! Tr)] (C8) 


Values of and are tabulated in table I for 7 equal 
to 0.074, 0.148, 0.222, and 0.296. The asymptotic magni- 
tudes of these functions are given by slender-body theory. 
Hence, one can easily show for large x 



Figure 43. — Definition of symbols used in equations (C3) and (C6.) 



742 


REPORT 1282 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Fiouse 44. — Definition of symbols used in equation (C7). 

Notice that both functions have essentially reached their 
asymptotic values for large x by the time 5= — 1 +5+0.5. 
By applying simple tabulative procedures to equation 
(C8) — for example, listing m £ and n { in reverse order and 
accumulating multiplications of adjacent terms — the value 
of (po? for any A 0 (x) representable by equation (C7) is readily 
calculated. 

The velocities induced by higher order multipoles can be 
calculated in a similar fashion. Because of the asymptotic 
behavior of the M’s and N’s, however, one is led into the 
numerically inefficient process of obtaining small numbers 
from differences of large numbers. For the velocities <p r and 
<pe, the following is a method for circumventing this difficulty. 

It follows from equations (9) or (19), that for small 
r, <p n - and <p„ e can be expressed in terms of the multipole 
strengths a n (x )— as defined by equation (7)— by the equa- 
tions 

r»+i / G 0 + Cir+. . .+Or’+. . . (ClOa) 
\cos r>BJ 

’~ H ■ ■+ D "'+ • ■ • <° 10b > 
where for v<n, n> 1 




\ - n u,v ( n 2) 

r K) 


(x), v even (Cl la) 


■ odd 


D v - 


"+s 

(— 1) 

A 4 X 



v even (Cl lb) 


v odd 


Consider now the velocities <p^ and <pi e induced by the 
multipole strength defined by equation (C3) ; thus 

-£^=mM ;t Jx)+iiN s £c) (Cl 2a) 


‘Pit 


where values of the M’s and N’s are listed in tables II and 
m. Their asymptotic values, as given by equations (CIO) 
and (Cll), are 




19+205 

1507JT 3 


N t 


1505 1 +2755+126 

1507IT 3 


M h 


19+205 
15 07+® 


150P+ 2755+ 126 1 | 

H 150 ifP 2,+J 


(Cl 3) 


and these are also given in the tables. 

As the tables show, equations (C13) are sufficiently 
accurate approximations to M and N for practical calcu- 
lating purposes when 5>— 1 +5+0.5. Hence, the velocities 
at the point 5,5 induced by the multipoles in the interval 
— 1<|<5— 5— 0.5 can be calculated using equations (C13). 
In terms of the distribution for + 2 (f) — which is equal to 
® 3 ®(f), see equation (16) — this means the multipole dis- 
tribution shown in the upper part of figure 45 can be cal- 
culated by means of the asymptotic formulas and the result 
added to that obtained for the distribution shown in the 
lower part of figure 45 by use of equations (Cl2b) and tables 
II and HI in a manner identical to the one represented by 
equations (C6), (C7), and (C8). 

The value of ^(5,5) induced by a multipole distribution 
such as that shown in the upper part of figure 45 is, on the 
basis of equations (CIO) and (Cll), 



(7+3+ 
7 rr 3 


2 



(Cl 4) 



r sin 28 


mM 3) (x)+nN le (p) 


(Cl 2b) 


FiGtrRE 45. — Range of application of equations (C12) and (C13). 
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whore p' t af (x)=Ai(x) and x t is shown in figure 45 and de- 
fined by the equation 


£<=0.4+2;— r 


[10(5— r)] 
10 


(C15) 


the symbol [10(5— r)] meaning, as before, the highest integer 
value contained in 10(5— 7). A similar result can be derived 
for ip-ij and one has finally for — 1 <[10(5— r)+ll]<5 


r sin 26 


cos 26 



and for [10(5— r) + ll]>6 


To i )+*- Ar ><( i -y)] ' 

(Cl6a) 


(Cl 6b) 


<P2g(x) 

[10(3g+H] 

r sin 28 

[lOff-^TH-llhS 


H(x,r) 


poet-n+ii] 

= 

cos 2 8~ 

[io(*4h-u}-6 


3(x,r)- 


(Cl7a) 

W i y) wt F( i Tr)] + 

(Pm 


The streamwise gradients of induced velocities can also be 
defined in terms of Mb and JV's as were the velocities 

<Pf and = Vg. Thus 

T 

‘po s =mH 0 2+nN 0i 


cos 29 


m'M^+nN^ 


Values of M 0 ~, N 0l , M s ,, and for r equal to 0.148, 0.222, 
and 0.296 are given in tables IV and V. 


APPENDIX D 

NUMERICAL METHOD USED TO CALCULATE BODY SHAPE 


The method used to calculate the body shape was a stand- 
ard step-by-step solution to the two simultaneous total 
differential equations (eq. (43) in the text) 


dx~V 0 R 1 


<1R 1 

dx~U„ 


(p T {x J R 1 Q') 


(Dl) 


The essentials of the process are recognized from the following 
computing-sheet heading set up for initial values of 9 and R 
equal to 30° and 0.148, respectively, where x=x/L 0 and 
R—PR/L 0 
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APPENDIX E 

OPTIMUM VARIATION OF a, (x) FOR A FIXED VALUE OF f° a, (x) dx 

J— L o 

If /(a :)=/(— x), then 


Di =~Ai£° L Jj [^.)+^i)i fe m (x a )+<(x 2 )] 


In 


Xi — Xj 


dxidxi 


(El) 


where the variation of a?' (x) is fixed, pose the restraint 

~ f a i (x)dx=F*(0)j constant (E2) 

U 0 J-L 0 4 

and ask for the function a 2 (a;) which minimizes D% for a given 
value of the constant. 


/(zi)/te)ln| a - ; — *2 1 dXidx 2 
/(a:i)/(a^)ln[a^— 

and since a^— A 0 ) = q^ u (— Z„) = a 2 <- ® (— _L o )=0 

f oc i {x)dx=—\ f x?a^> (x)dx 
J-l 0 6 J-l c 

Therefore, the standard variational problem 


(E3) 
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reduces to 

5 { 2 XUS f-L f-L ^ (Xt) + “^ (X ‘ )] ^ +0f> (XS)] 

^ ^ 1 chidx.-^j- J ^ xfe® (x 1 )dx 1 J- =6 


In. 

and this becomes 


LI 


£. ‘ [^] *■ {m, £. wi 


In. 


xl—xl 


Li 


dx* 


“T } =< 


Integrating three times by parts, using the relations 
5^(~L 0 )=5c4» (-i 0 )=5a® (-£„)= 0 
5a 2 (0)=5a 2 ® (0)=0, by symmetry 

and 

'0 


lim ~ f° [«,*(**)+<(*,)] 

n -+ 0 OX I J -L' J^o 


dx 2 


-,i»2x, f glaM to) fe=0 
1,-0 J-i * 1— *2 


yields 


[dxj+C/^A,^ J" =0 


In. 


xl—xl 


Li 


where 


02 


(E4) 


By the fundamental lemma of the calculus of variations 
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One can also show 
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dx,— TJ^x{ J- =0, 

Lfy^>Xi~^>0 


Hence, 

0 «£ (x-j) 


j: 


Uo 


In 


xl—xl 


Li 


dx i =Y () +72xf — Kzf 
'° <h m (x i ) 


xf— xf 



J-x„ u 0 H 

1 Li 

Integrating by parts and changing the notation 

x\ 

Vi=jV 



X? 

Vl ~Li’ 




dx 2 (E5) 


p/w^ =T<+7wi _ v ^_ f-sM 

Jo ni — vt Jo vi 


1 g(vt)dv2 

Vi 


(E6) 


Equation (E6) is the familiar singular integral equation 
known, in aerodynamic applications, as the airfoil equation. 
Its inversion is discussed, for example, in reference 3. If one 
solves equation (E6) and applies the conditions 


then 


7o= 


o®(0)=o 2 ®(L o )=0 


4 \>L* 48 \ 0 L 0 

’ Ti== -16£- 


15w 


and 


^-^{[K0-]VK0 


+ 


of X V It, L—^Li — x'X 

\L.) 1 ’ > i.+V3YV>/ 

(E7) 


(E8) 


Now set 


so that ^*(0) is a known constant. Then 


2 X„L» 


^-Lrvr =F ^ 0) l Ri — F ** (p) (E9) 

Using the above expressions, one can show 


ajjx) , a*(x) 15 ttR* c x 


U 0 1 U„ 


32 Li 


[^*(0)-W(0)] 


{K0-#-(0+K0-^=:> 

(E10) 

The wave drag can be calculated by combining equations 
(E5) and (El). 


\dx - ! 


Integrate three times by parts and there results, finally 

D * = W qL>3 (tJ [E*(0)+I£W (Ell) 


REFERENCES 

1. Heaslet, Max. A., Lomax, Harvard, and Spreiter, John R.: Lin- 

earized Compressible-Flow Theory for Sonic Flight Speeds. 
NACA Rep. 956, 1950. 

2. Lamb, Horace: Hydrodynamics. Sixth ed., Dover Publication 

(New York), 1945 

3. Lomax, Harvard, Heaslet, Max. A., and Fuller, Franklyn B.: 

Integrals and Integral Equations in Linearized Wing Theory. 
NACA Rep. 1054, 1951. 

4. Watson, G. N.: A Treatise on the Theory of Bessel Functions. 

Second ed., Cambridge Univ. Press (Cambridge, England), 1962. 



THE REDUCTION OF WAVE DRAG OF WING AND BODY COMBINATIONS 


745 


5. Erd61yi, Arthur: Table of Integral Transforms. VoL I, McGraw- 
Hill Book Co., Inc., 1954. 

0. Hayes, Wallace D.: Linearized Supersonic Flow. North American 
Aviation, Inc., Rep. No. AL-222, June 1947. 

7, Lomax, Harvard, and Heaslet, Max. A.: Recent Developments in 

the Theory of Wing-Body Wave Drag. IAS Preprint 617, 1956. 

8. von KArmdn, Th.: The Problem of Resistance in Compressible 

Fluids. (Fifth Volta Congress) Roma Reale Accademia D’ltalia 
1936. 

0. Webster, Arthur Gordon: Partial Differential Equations of Mathe- 
matical Physics. Second ed., Hafner Publishing Co., Inc. (New 
York), 1950. 

10. Niolsen, Jack N., and Pitts, William C.: General Theory of Wave- 
Drag Reduction for Combinations Employing Quasi-Cylindrical 
Bodies With an Application to Swept Wing and Body Combina- 
tions. NACA TN 3722, 1956. (Formerly NACA RM A55B07). 


11. Jones, Robert T.: Theoretical Determination of the Minimum 

Drag of Airfoils at Supersonic Speeds. Jour. Aero. Sic., vol. 19, 
no. 12, Deo. 1952, pp. 813-822. 

12. Jones, Robert T.: Theory of Wing-Body Drag at Supersonic 

Speeds. NACA RM A53H18a, 1953. 

13. Lomax, Harvard: The Wave Drag of Arbitrary Configurations in 

Linearized Flow as Determined by Areas and Forces in Oblique 
Planes. NACA RM A55A18, 1955. 

14. Nielsen, Jack N., and Pitts, William C.: Wing-Body Interference 

at Supersonic Speeds With an Application to Combinations With 
Rectangular Wings. NACA TN 2677, 1952. 

15. Jones, Robert T.: Some Recent Developments in the Aerodynam- 

ics of Wings for High Speeds. Zeitschrift Ftir Flugwissenschaf- 
ten, 4. Jahr., Heft 8, Aug. 1956, pp. 257-262. 

16. Graham, Ernest W.: The Pressure on a Slender Body of Non- 

Uniform Cross-Sectional Shape in Axial Supersonic Flow. 
Douglas Aircraft Co., Inc., Rep. No. SM-13346-A, July 1949. 



746 


REPORT 1282 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TABLE I.— VALUES OF M o7 AND N a7 . 


TABLE II— VALUES OF M r , ; AND Nxf. 



TABLE m— VALUES OF AND N, r 


TABLE IV— VALUE OF M 0 -, AND N„. 




Exact Asymp- Exact Asymp- 
totlo totio 


0 

6.62 

13.69 
18.38 
23.48 

23.70 
33.92 
39.16 
44.47 I 44.41 


49.63 49.64 

64.86 64.88 





2.22 7.69 7.72 

2.41 9.33 9.35 

2. 61 1L 11 1L 14 


Exact Asymp 
totio 


-a 852 

0 

-.802 

L 86 

752 

2.93 

-. 702 

3.34 

-.652 

3.93 

-.602 

4.58 

—.552 

5.22 

-.502 

5.89 

-.452 

6.49 

-.402 

7.16 

-.352 

7.83 



-a 852 

0 

0 

-.802 

- 1.663 

-.442 

-.752 

2.529 

- 3.024 

-.702 

.369 

- 1.010 

-.652 

.191 

-.723 

-.602 

.121 

-.572 

-.552 

.085 

-.477 

-.502 

.063 

-.410 

-.462 

.049 

-.361 

-.402 

.040 

-.323 

-.352 

.033 

-.292 

-.302 

.027 

-.267 

-.252 

.023 

-.245 

-.202 

.020 

-.227 

-.152 

.018 

-.212 

-. 102 

.015 

-.199 

-.052 

.014 

-.187 

-.002 

.012 

-.176 

. 048 

.011 

-. 167 

.098 

.010 

-.168 

.148 

.009 

-.151 

.198 

.008 

-. 144 

.243 

.007 

-. 138 

.298 

.007 

-.132 

.343 

.006 

-.127 

.393 

.006 

-.122 

.448 

.005 

-.117 




-.578 

.165 

-.020 

-.504 

-.523 

.099 

-.498 

-.454 

-.478 

.070 

-.420 

-.404 

-.423 

.053 

-.365 

—.354 

-.378 

.042 

-.324 

-.304 

-.323 

.034 

-.202 

-.254 

-.273 

.023 

-.266 

-.204 

-.223 

.024 

-.243 

-.164 

-. 178 

.020 

-.220 

-. m 

-.123 

.018 

-. 210 

-.054 

-.078 

.016 

-. 107 

-.004 

-.023 

.014 

-. 185 

.040 

.022 

.012 

-. 174 

.090 

.072 

.011 

-.105 

.146 

.122 

.010 

-. 167 

.100 

.172 

.009 

-.149 

.240 

.222 

.003 

-. 142 

.296 

.272 

.007 

-.130 

.340 

.322 

.007 

-.131 

.390 

.372 

.006 

-. 125 

.440 

.423 

.006 

-.120 






4.82 4.85 

6.67 6.69 


































































































































































